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Zusammenfassung. Die kosmologische Information, die im Signal des schwa- 
che-Gravitationslinsen-EfFekts verborgen ist, lasst sich mit Hilfe des Potenzspek- 
trums der sogenannten Konvergenz analysieren. Wir verwenden den Fisher-In- 
formation-Formalismus mit dem Konvergenz-Potenzspektrum als Observable, um 
abzuschatzen, wie gut zukiinftige Vermessungen schwacher Gravitationslinsen die 
Expansionrate und die Wachstumsfunktion als Funktionen der Rotverschiebung 
einschranken konnen, ohne eine bestimmtes Modell zu deren Parametrisierung zu 
Grunde zu legen. Hierfiir unterteilen wir den Rotverschiebungsraum in Bins ein 
und interpolieren diese beiden Funktionen linear zwischen den Mittelpunkten der 
Rotverschiebungsbins als Stiitzpunkte, wobei wir fur deren Berechnung ein Be- 
zugsmodcll verwenden. Gleichzeitig verwenden wir diese Bins fur Potenzspektrum- 
Tomographie, wo wir nicht nur das Potenzspektrum in jedem Bin sondern auch 
deren Kreuzkorrelation analysieren, um die extrahierte Information zu maximie- 
ren. 

Wir stellen fest, dass eine kleine Anzahl von Bins bei der gegebenen Prazision 
der photometrischen Rotverschiebungsmessungen ausreicht, um den Grofiteil der 
Information zu erlangen. AuBerdem stent sich heraus, dass die prognostizierten 
Einschrankungen vergleichbar sind mit derzeitigen Einschrankungen durch Beob- 
achtungen von Clustern im Rontgenbereich. Die Daten von schwachen Gravitati- 
onslinsen alleine konnten manche modifizierte Gravitationstheorien nur auf dem 
2a Niveau ausschliefien, aber wenn man den Prior von Vermessungen der kosmi- 
schen Hintergrundstrahlung beriicksichtigt, konnte sich dies zu einem 3a Niveau 
verbessern. 

Abstract. The cosmological information encapsulated within a weak lensing 
signal can be accessed via the power spectrum of the so called convergence. We 
use the Fisher information matrix formalism with the convergence power spectrum 
as the observable to predict how future weak lensing surveys will constrain the 
expansion rate and the growth function as functions of redshift without using 
any specific model to parameterize these two quantities. To do this, we divide 
redshift space into bins and linearly interpolate the functions with the centers of 
the redshift bins as sampling points, using a fiducial set of parameters. At the 
same time, we use these redshift bins for power spectrum tomography, where we 
analyze not only the power spectrum in each bin but also their cross-correlation 
in order to maximize the extracted information. 

We find that a small number of bins with the given photometric redshift mea- 
surement precision is sufficient to access most of the information content and that 
the projected constraints are comparable to current constraints from X-ray cluster 
growth data. This way, the weak lensing data alone might be able to rule out 
some modified gravity theories only at the 2a level, but when including priors 
from surveys of the cosmic microwave background radiation this would improve to 
a 3a level. 
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Chapter 1 
Introduction 



1.1 Motivation 

It is a truth universally acknowledged that any enlightened human being in 
possession of a healthy curiosity must be in want of an understanding of the 
world around him. In the past few decades, two mysterious phenomena have 
emerged in cosmology that challenge our understanding of the Universe in 
the post provoking way: dark matter and dark energy. The "stuff" that we 
observe, that makes up the galaxies and stars, the protons and the neutrons, 
the planets and the moons, the plants and the animals and even ourselves, 
everything we ever knew and all of what we thought there ever was, makes up 
only a meager 4% of the cosmos. At first glance, this presents a huge problem 
for cosmology. But in the grand picture, it might be an opportunity for all 
of physics to make progress. Now, for the first time in history, cosmology 
is inspiring particle physics via dark matter for the hunt of a new particle, 
bringing physics to a full circle. At the same time, dark energy is causing 
embarrassment in the quantum field theory department by deviating from 
the predicted value by 120 orders of magnitude. Closing all gaps in our 
knowledge about the Universe on a large scale is therefore necessary in order 
to understand physics and our world as a whole. 

Experiments are the physicist's way of unlocking the Universe's secrets, 
but they are incredibly limited in the field of cosmology. The studied object is 
unique and cannot be manipulated, only be observer through large telescopes 
either on the ground or in space. There is no shortage of theories that could 
explain some of the issues that are troubling cosmologists these days, but 
they all look very similar and it is hard to tell them apart by looking at the 
sky. Subtle clues need to be picked up to succeed, and for this we need better 
and better surveys. To shed light onto the dark Universe, several missions 
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are being planned, designed, or are already in progress. In the design stage 
in particular, it is interesting to see what kind of results one can expect 
depending on the experiment's specifications. 

A lot of dark energy models are degenerate in the sense that some of 
their predictions match or are indistinguishable by reasonable efforts. The 
most popular ones predict critical values for some parameters, so unless we 
are provided with a method that allows for infinite precision in our mea- 
surements, we will never know, for instance, whether the Universe is flat 
or has a tiny, but non-zero curvature that is too small to be detectable, or 
whether dark energy is really a cosmological constant or is time-dependent 
in an unnoticeable way. We hope to break at least some of the degeneracies 
by investigating how growth of structure in the matter distribution evolved 
throughout the history of the Universe. 

Information about how the matter used to be distributed is naturally ob- 
tained by analyzing far away objects, since the finite speed of light allows 
us to look as far into the past as the Universe has been transparent. But 
the nature of dark matter is hindering us from getting a full picture of the 
matter distribution, because only 20% of all matter interacts electromagnet- 
ically. The only sensible way to detect all of matter is via its influence on 
the geometry of space-time itself. Distortions of background galaxies will 
reveal perturbations in the fabric of the cosmos through which their light 
has passed, unbiased in regard of matter that happens to be luminous or 
not. In particular, the signal coming from the so-called weak lensing effect 
is sensitive to the evolution of both dark energy and growth related parame- 
ters. To make the most out of this phenomenon, we "slice up" the Universe 
into shells in redshift space and study the correlation of these redshift bins, 
a technique which has been dubbed weak lensing tomography. 1 

1.2 Objective 

By assuming the validity of Einstein's general relativity, all observations are 
biased towards the standard model and the growth of structure might not 
have been measured correctly. Instead, we want to allow for the possibility of 
modified gravity. Our goal is to investigate how future weak lensing surveys 
like Euclid can constrain the expansion rate H(z) and the growth function 
G(z) without assuming a particular model parameterizing those two quan- 
tities. We first select a suitable number of redshift bins in which we divide 
a given galaxy distribution function. Starting from a fiducial model with a 
list of cosmological parameters that take on particular values determined by 



1 Tomos (rd^o<j): Greek for slice, piece. 
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Wilkinson Microwave Anisotropy Probe (WMAP) 2 , we then take the values 
of H and G in each bin as a constant, treat them as additional cosmological 
parameters, and rebuild these two functions as linear interpolations between 
supporting points in the redshift bins. 

Using these new functions, we calculate the weak lensing power spectrum 
in each bin as well as their cross-correlation spectra based on the non-linear 
matter power spectrum, which we take from fitting formulae found by other 
groups. The weak lensing power spectrum is then used in the powerful Fisher 
matrix formalism, which allows us to estimate the uncertainty of all cosmo- 
logical parameters, including the values of H and G at the chosen redshifts. 
With these forecast error bars, we can compare competing theories of modi- 
fied gravity to our simulated results of next generation weak lensing surveys. 

1.3 Symbols and notation 

We use units where the speed of light equals unity. Vectors are lower case and 
printed in bold (e.g. x), while matrices are upper case and printed sans serifs 
(e.g. F). The metric has signature — h++. Derivatives are sometimes written 
using the comma convention: dcfi/dx 1 = <p^. The logarithm to base ten is 
denoted by lg, and the logarithm to base e by In. The Fourier transform 
of a function f(x) is written with a tilde: f{k). Due to the one-to-one 
correspondence of the redshift and the scale factor, a = 1/(1 + 2;), we may 
denote the dependence on either one of those quantities interchangeably, for 
instance it is obvious that H(a) = H(l/(l + z)), so we might as well consider 
H(a) as a function of z and write H(z). Sometimes, the value of redshift 
dependent quantities like fl m (z) as of today is given as just fl m , for instance. 
A list of commonly used symbols follows. 



a 


Scale factor 


X 


Complex ellipticity 


Cij 


Covariance matrix 


Sij 


Kronecker delta 


Sb 


Dirac delta function 


S m (z) 


Matter density contrast 


A 2 


Photometric redshift error 


E(z) 


Dimensionless Hubble parameter, E(z) = H(z)/H 


F 


Fisher matrix 


f(z) 


Growth rate 



2 A probe of the cosmic microwave background radiation, see 
|http : / /ma p . gsf c . nas a . gov[ 
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G Newton's gravitational constant 

G(z) Growth function 

7 Growth index 

7i n t Intrinsic galaxy ellipticity 

7i Shear, i — 1,2 

7 PP n Parameterized post-Newtonian parameter 

hi,gi Values of the logarithm of H(z),G(z) at redshift z i: i = 

Hq Hubble constant 

h Dimensionless Hubble constant, h = iJn/100km/s/Mpc 

Jij Jacobian matrix 

k Wave number 

k Convergence 

i Multipole 

C Likelihood, Lagrangian 

n(z) Galaxy distribution function 

Tii(z) Galaxy distribution function for the i-ih redshift bin 

rii Average galaxy density in the i-th redshift bin 

n s Scalar perturbation spectral index 

ritf Angular galaxy density 

M Number of redshift bins 

u x Reduced fractional density for component x: u x = Q x h 2 

fi x (a) Fractional density for component function of the scale 

factor or redshift (x = m,b,c for matter, baryons, cold dark 

matter) 

fi m Fractional matter density today 

P m (k) Matter power spectrum 

P(£) Convergence power spectrum 

p Pressure 

P(A\B) Probability of A given B 

$, ^ Scalar gravitational potentials 

r(z) Comoving distance 

p Density 

o"8 Fluctuation amplitude at 8Mpc/h 

cr(x) Uncertainty of the quantity x 

t Reionization optical depth 

t Cosmic time 

6 Vector in parameter space; angular position 

w Equation-of-state ratio 

W(z) Window function 

z Redshift 
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z m Median redshift 

Zi, h Endpoint, center of the z-th redshift bin 



Chapter 2 

Theoretical preliminaries 



This section outlines the basic theory behind dark energy by giving a brief 
historical introduction. Furthermore, the basics of weak lensing and the 
Fisher information matrix formalism are explained, as well as the concept 
of power spectra, an immensely important tool not only in the remainder 
of this thesis, but in all of modern cosmology. A much more detailed treat- 
ment of the topics presented in this section can be found in many excellent 
standard text books, such as Carroll (2003) for the derivation of cosmol- 



ogy from general relativity, Weinberg (2008) and Amendola and Tsujikawa 



(2010) (hereafter referred to as A&T) for an up to date reference to mod- 
ern cosmology, dark energy and useful statistical tools, Peebles (1980) also 
for the statistics and Bartelmann and Schneider (1999) for a comprehensive 
treatment of weak gravitational lensing. 



2.1 Dark energy: A historical summary 

A short while after Albert Einstein derived his famous field equations he 
proposed the idea of inserting a cosmological constant originally in order to 
allow for a non-expanding Universe, because he felt that a world that was 



spatially closed and static was the only logical possibility (Einstein, 1917). 
This view was the general consensus at the time, considering that all visible 
stars seemed static and other galaxies had not been discovered yet. From 
the first and second Friedmann equations 



d\ 2 8tiG k 



a J 3 a 2 
a __ A-kG 
a 3 



P-~ 2 , (2-1) 
(p + 3p) , (2.2) 
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we can easily see, that for a static solution with a = a 
component in the Universe that bears a negative pressure 



0, we need a 



V 



with an overall positive spatial curvature that is finely tuned to 



k 8ttG 



-P- 



(2.3) 



(2.4) 



a? 3 

At this point, we mention that it is often useful to plug the time derivative 



of eq. (2.1) into eq. (2.2) to arrive at the continuity equation 



p=-2-(p + p). 
a 



(2.5) 



It can be shown that for dust, i.e. non-relativistic baryonic or dark matter, 
the pressure vanishes and for radiation or ultra-relativistic matter the pres- 
sure is one third of the energy density. This means that ordinary forms of 
energy cannot produce the desired result. Einstein recognized that we the 
Lagrangian in the Hilbert-Action can be changed to 



S 



1 



16nG 



^x^g{R - 2A) 



(2.6) 



by introucing a constant A. We can interpret this term as a form of energy 
density p\ = A/8nG that does not depend on the scale factor, which leads to 
p = —p according to eq. (2.5). Thus p in eqs. (2.1) and (2.2) can be replaced 
by Pm + Pa- Then a Universe in which p m = 2p\ or A = 47iGp m satisfies eq. 



(2.3). 



However, Einstein's conservative view of the world did not live up to 
experimental facts. Only a short time later, after Edwin Hubble's revolu- 
tionary discovery of apparently receding galaxies in 1923, Einstein readily 
discarded the idea of a static Universe. Because an expanding Universe does 
not necessarily need a cosmological constant (see fig. 2.1), the introduction 



of the cosmological constant was considered a mistake for decades after that 



(Peebles and Ratra, 2003) 



It was not until the early 1960s, when A had to be revived to explain 
a new measurement of H by Allen Sandage that was more accurate than 
the original one by Hubble by one order of magnitude. After a proper re- 
calibration of the distance measure, he found Hq w 75km/s/Mpc, causing 
an incompatibility of what was then concluded to be the age of the Uni- 



18 



Chapter 2. Theoretical preliminaries 



verse (7.42 Gyr) with that of the age of the oldest known star in the Milky 



Way, which was at the time estimated to be greater than 15 Gyr (Sandage 



1961). Even though the numbers were slightly off according to our current 



knowledge, Sandage was on the right track and immediately suggested that 
a cosmological constant could easily resolve this issue. This was the first hint 
that the Universe was not only expanding, but expanding in an accelerated 
manner. 

Several decades later, the highest authority in physics — the experiment- 



put an end to all speculations. In the late 1990s Riess et al. ( 1998 ) observed 



ten supernovae of type la that they used as standard candles to infer their 
luminosity distance as a function of redshift and found that they were on 
average 10%-15% farther away than one would expected in an open, mass 
dominated Universe. They were able to rule out the = case at the 9a 
confidence level, leading the way for precision cosmology. Almost simulta- 



neously, an independent measurement by Perlmutter et al. (1998) confirmed 



their findings by analyzing 42 high-redshift supernovae of type la. Thus, the 
notion that the expansion of the Universe was accelerating was confirmed. 
Our failure to understand the nature of this accelerated expansion is re- 
flected in the name for this phenomenon: dark energy. It is important to 
stress that the underlying physical cause of the acceleration does not have 
to be a new form of energy, it might as well be a new form of physics that 
mimics the effects of a cosmological constant, which is often misunderstood 
even by physicists unfamiliar with cosmology. Whether we are faced with 
a "physical darkness or dark physics" , to say it with the words of Huterer 
and Linder (2007), is the whole crux of the matter. Evidence independent 



of the distance-luminosity relation for type la supernovae which supports 
the idea of acceleration has steadily increased since the findings of Riess and 



Perlmutter. For instance, the late integrated Sachs- Wolfe effect (Scranton 



et al, 2003), or the cosmic microwave background radiation (Komatsu et al. 



2010) and a dark energy component is part of today's standard model of the 



Universe. But does this mean that cosmology as a science is done? 

On the contrary. Even though the so called flat ACDM model (cold 
dark matter with a cosmological constant) with only six free parameters 
(u m ,Ub,QA,n s ,<J8,T) 1 agrees remarkably well with all observations, most physi- 
cists are still not satisfied with it for a number of reasons. First of all, the 
value of the cosmological constant is incredibly small. When we express 
A in natural units, i.e. where the Planck length lp = y/hG = 1, we get 
A = 3.5 x 10~ 122 . Even worse, when we interpret the cosmological constant 
simply as the vacuum energy predicted from quantum field theory, the value 



1 The exact meaning of those parameters will be explained in section 3.4 
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-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 
f 



Figure 2.1 Solutions of the Friedmann equation demonstrating 
the evolution of the Universe. We have set the Hubble constant Hq 
to unity for convenience, which implies that t = — 1 corresponds 
to the Hubble time — 13.9Gyr. All curves pass t = (today) 
with a(0) = d(0) = 1. Top panel: From the top: The expansion 
history and future for our Universe with (O m = 0.3, Q\ = 0.7), 
an open matter-only Universe (0.3, 0) and the Einstein-de Sitter 
model (1,0). Bottom panel: Some more exotic solutions: A pure 
cosmological constant (0,1), the Milne model (0,0) and a closed 
Universe (6, 0). 
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diverges. Because it is widely believed that new physics are needed at the 
Planck scale where neither quantum effects nor gravity are negligible, namely 



a theory of quantum gravity (Carroll et at, 1992), one can renormalize the 



integral over all modes of a scalar, massive field (e.g. a spinless boson) with 
zero point energy ftw/2 by introducing a cutoff at Planck length. This gen- 
erates a finite value, but one that is still 120 orders of magnitude too large, 
making it arguably the worst prediction in the history of physics. Finding a 
mechanism that would cause fields to cancel out each other's contributions 
to the vacuum energy density would be a challenge, but a mechanism that 
leaves a tiny, positive value seems completely out of reach for now. The 
most promising resolution might lie in super symmetry, super gravity or su- 
per string theory. This is commonly referred to as the cosmological constant 



problem (Weinberg, 1989). 



The second problem is the coincidence problem, which describes the cu- 
rious fact that we live exactly in a comparably short era in which matter 
surrenders its dominance to dark energy. This can be visualized if we con- 
sider the time dependence of the dimensionless dark energy density 2 : 



A (a) 



Pa 



Pa 



1 



1 



Pcrit(a) p m (a) + p A 



PA 



(2.7) 



because p m scales as a 3 while p\ stays constant, when an is the scale factor 
today, usually normalized to unity. Plotting this with its derivative on a 



logarithmic scale (see fig. 2.2) shows the coincidence. 

Until now, there is no satisfactory solution to these problems. String the- 
orists argue that we might live in one of 10 500 realizations of the Universe, 
corresponding to the number of false vacua that are allowed by string the- 



ory (Douglas, 2003). This approach comes down to an anthropic argument, 



stating that it is no surprise that the cosmological constant has the value 
that it has, since any other value might be realized in a different Universe, 
but would not lead to the growth of structure required for sentient life that 



can ponder the value of natural constants (Susskind, 2005). This argument 



is controversial amongst physicists. As one of the leading cosmologists Paul 



Steinhard puts it in a Nature article (Brumfiel, 2007): "Anthropics and 
randomness don't explain anything." 

Alternative theories include the postulation of Quintessence, a new kind 
of energy in the form of a scalar field with tracker solutions mimicking a 



cosmological constant (Zlatev et al, 1999), or a modification of Einstein's 



2 Radiation energy scales as a 4 and was only relevant at very early stages of the 
Universe, which is why we shall not consider it here. 
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Figure 2.2 Visualization of the coincidence problem. Shown here 
is the dimensionless dark energy density £l\(a) in dependence of 
the scale factor a and its derivative on a logarithmic scale. High- 
lighted with dashed lines are the times of big bang nucleosynthesis 
(BBN) at a ~ 10 -10 , recombination at a ~ 1/1000 and today at 
a = 1. Dark energy started to dominate right when structures 
emerged in the Universe. 



theory of general relativity, where, in the most popular case, the Ricci scalar 



in general relativity is replaced as the Lagrangian by some function f(R) (De 



Felice and Tsujikawa, 2010). Unfortunately, current observational constraints 



are insufficient when it comes to discriminating between competing theories. 



2.2 Weak lensing 

When light passes through a gravitational potential, its trajectory is bent in a 
way described by general relativity. This process, dubbed "lensing" , can alter 
the shape, size and brightness of the original image. We usually distinguish 
between strong, weak and micro lensing. The mechanism of micro lensing 
relies on small objects of low magnitude being the lens, such as brown or white 
dwarfs, neutron stars, planets and so on, that transit a bright source and 
temporarily increase the source's brightness on time scales of several seconds 



to several years (Paczynski, 1996). The other two cases are static on those 
time scales, since more massive lenses like galaxies or clusters are involved, 
which means that the distances involved are several orders of magnitude 
larger. Strong lensing occurs when multiple images of the same galaxy from 
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behind a super massive object can be observed, as seen in the bottom left 
corner of the left panel in fig. 2.3| 



>■ .■• •■ «! a V' ^r . yb ■ ■- E 




Figure 2.3 ie/i panel: Simulation of a gravitational lens with a 
case of strong lensing in the lower left corner at an average redshift 
of one. Right panel: Enhancement of the upper right section of the 
left panel, depicting lensing in the weak regime. The contours are 
ellipses derived from the quadrupole moment of each galaxy image. 
The two bars in the upper right corner show the average shear, 
where the lower one represents the expected shear and the upper 
one the actual shear. The deviation stems from shot noise due 
to the intrinsic ellipticity and random orientation of the galaxies. 



Taken from Mellier (1998) 



While strong and micro lensing are relatively rare, weak lensing is an ef- 
fect that can be observed over large areas of the sky (see upper right corner of 
the right panel in fig. 2.3). Its power lies in statistical analysis, since weakly 



lensed objects are only slightly distorted and impossible to detect individu- 
ally. Only the average distortion field reveals the cosmological information 
that we are looking for. As it turns out, weak lensing changes the ellipticity 
of galaxies (in a first order approximation), but the intrinsic ellipticity al- 
ways dominates the shape. We need to gather large enough samples and then 
subtract the noise, which is relatively simple if the magnitude and direction 
of the intrinsic ellipticity is uncorrelated to the signal. This is unfortunately 
not entirely the case due to tidal effects (intrinsic alignment), and this is 
only one of the many challenges that weak lensing harbors. The most obvi- 
ous one is illustrated in fig. 2.4 Others include photometric redshift errors, 



calibration errors and uncertainties in power spectrum theory. A lot of these 
systematic errors can be accounted for by introducing several nuisance pa- 



2.2 Weak lensing 
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rameters (Bernstein, 2009). The trade-off is that a high number of nuisance 
parameters diminish the merit of the Fisher matrix formalism, as there are 
degeneracies to be expected, so we will only account for photometric redshift 
error and the intrinsic ellipticity in this thesis. 

a) b) c) d) e) 

Figure 2.4 A schematic illustration of the process of weak lens- 
ing detection, demonstrating its difficulty, a) The original source 
galaxy, b) Weak lensing distorts the shape and adds a shear to the 
image, c) The image is convolved by the telescopes point spread 
function and, in case of ground based surveys, the atmosphere, d) 
The effect of the finite resolution of the detector, e) Additional 
noise is applied to the image. The challenge now is to infer b) 
from e). (Image of M31 by courtesy of Robert Gendler.) 

An expression describing the shear can be obtained by perturbing the 
Friedmann-Lemaitre-Robertson- Walker metric in the Newtonian gauge 

ds 2 = -(1 + 2^)dt 2 + a 2 (t)(l + 2$)(dx 2 + dxj + dx 2 3 ) , (2.8) 

where \I/ and $ are scalar gravitational potentials. We can then solve the 
geodesic equation for a light ray 

— +^ = 0, (2.9) 

which can be rewritten in our case as 

= (2-10) 

Thus, we obtain the lensing equation 

e i = Oi + J ^'^1-^^(^01,^02,0 (2.11) 

with the lensing potential ip = $ — Here r is the radial comoving coordi- 
nate, 9i = Xi/r is the angle of the light ray with respect to the r-axis, and 
(xi,x%) are displacement coordinates perpendicular to the r-axis. 

The distortion of a source image at distance r s to first order is a linear 
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Observer 




Source plane 



Figure 2.5 Description of the lensing configuration. Note that 
comoving distances cannot be added trivially, i.e. r s ^ r j + n s . 



transformation described by the symmetric matrix 

d9„ 



A 



00, 



Oj 



with 



dr' I 1 ) r'^, 



Here we introduced the convergence 

1 



-«wl - 7i ~72 
-72 -«wl + 7i 



dr' 1-- r'(V,n+^22) 



(2.12) 



(2.13) 



which is a measure for the magnification of the image, and the shear 

1 rr ' 



(2.14) 



7i = -; 



72 = - 



dr' 1 r'(-0,n - ^,22) 



dr' 1 



12 



(2.15) 
(2.16) 



which describes the distortion. These quantities will become important when 
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we want to describe the cosmic shear by statistical means, in particular its 
power spectrum. 

We can measure the ellipticity of real astronomical images by computing 
their quadrupole moment, which is defined as 

Oij = J 1(6)9^ d 2 6, (2.17) 

where 1(0) is the luminous intensity of the galaxy image with its center at 
= 0. The complex ellipticity is then given by 

X= ; • 2.18 

5ll + 522 

It is straight forward to show that according to this definition the complex 
ellipticity for an ellipse with semiaxes a and b, rotated by an angle 0, is 

n 2 - h 2 
cr + tr 

To first order, a weak lens distorts a spherical object into an ellipse with a 
simple relation between the shear and the complex ellipticity, namely 

x = 2( ll+ i l2 ). (2.20) 

We can use this to compute the power spectra P 7i (£), which describe the auto- 
correlation of the shear field at the multipole £, the adjoint variable to 6. It 
can be shown that the convergence power spectrum P is related to the matter 
power spectrum and can be written to first order as a linear combination 
P = Pe = ciP 7l + C2-P 72 . Another linear combination, Pb = ciP 7l — c 2 P 72 , 
must vanish. These are usually referred to as the electric and magnetic part 
of the shear field. Thus, the convergence power spectrum contains all the 
information while the magnetic part provides a good check for systematics. 

But, as mentioned before, not all galaxies appear circular even if their 
intrinsic shape was a perfectly circular disc and there was no gravitational 
lensing distorting our view of the sky, as they are randomly oriented and thus, 
viewed from the side, look like ellipses. This intrinsic ellipticity, denoted by 
7i n t, can be reflected in a noise term that is added to the power spectrum, if 
we assume that the noise is uncorrelated to the weak lensing signal 3 : 

P = Psig„al + 7int^ 1 ) (2.21) 



3 As mentioned before, because of tidal effects this is not necessarily the case. 
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where n is the average galaxy density (A&T, ch. 14.4). 

The shear field is a vector field of dimension two. To access the full 
power of weak lensing, we must also include the third dimension, which can 
be done by considering the redshift of each source galaxy. Weak lensing 
works linearly, so we can add up the transformation matrices of all galaxies 
and find for the full transformation matrix 



dr'n(r')Dij(r') 



(2.22) 



where n(r) dr is the number of galaxies in a shell with width dr and the 
galaxy distribution function n(r) is normalized to unity. Making use of the 
fact that for any real, integrable functions f(x) and g{x,y) the identity 



dxf(x) / dyg{x,y) 



dy / dxf(x)g(x,y) 



(2.23) 



holds, we can rewrite eq. (2.22) to get 



Vij = / drw(r)ip :ij . 
Jo 

Here we abbreviated the weight function as 

(/■( / ) = / dr 1 1 — 



(2.24) 



rn(r 



(2.25) 



Thus, the total convergence field from eq. (2.14) now reads 



drw(r) (^n + tp. 



22) 



(2.26) 



By grouping all observed galaxies into redshift bins, we are able to extract 
more information out of the given data by calculating the power spectrum 
in each each redshift bin as well as their cross correlation. This is called 
power spectrum tomography (Hu, 1999). We shall cover this in more detail 
in section EH1 



2.3 Power spectra 

Almost every book and every lecture on modern cosmology begins with the 
cosmological principle: the assumption that the Universe is, simply put, 
homogeneous and isotropic. Now obviously the average density, say, 100 m 
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below the University Square in Heidelberg is quite different than the average 
density 100 m above, and it is still quite different at the midpoint between 
the Sun and Alpha Centauri - we say that the Universe has structure. But 
if we look at the Universe at sufficiently large scales of > 1 Gpc, the overall 



success of the standard model confirms our initial assumption (Komatsu 



et al, 2010). While there are research groups who investigate the possibility 



of genuine large scale inhomogeneities (Barausse et al. 2005 Kolb et al 



2005), we remain confident that there is good reason to believe that the 
cosmological principle applies to the Universe as a whole. 

The matter distribution in tour Universe consists mostly of structures 
at different scales, which can be roughly classified as voids, super clusters, 
clusters, groups, galaxies, solar systems, stars, planets, moons, asteroids and 
so on. This "lumpiness" is inherently random, but like all randomness it can 
be characterized by an underlying set of well defined rules by using statis- 
tics. The matter distribution in our Universe is one sample drawn from an 
imaginary ensemble of Universes with an according probability distribution 
function (PDF), and we are now facing the challenge of determining that 
PDF. 

We require our models to make precise predictions regarding the structure 
of matter. Naturally, no model will be able to predict the distance between 
the Milky Way and Andromeda for instance, but it should be able to predict 
how likely it is for two objects of that size to have the separation that they 
have. To quantify and measure the statistics of the irregularities in the matter 
distribution in the cosmos, we define the two-point correlation function £ as 
the average of the relative excess number of pairs found at a given distance. 
That is, if n,-^ i e {1,2}, denotes the number of point-like objects (nucleons, 
or galaxies, or galaxy clusters) found in some volume element dVi at position 
rj with a sample average (n,) = n dVi, we can write the average number of 
pairs found in dVi and dV 2 as 

M = n 2 Q dV x dV 2 [l + r 2 )}, (2.27) 

where no = N/V is the mean number density. If the matter distribution was 
truly random, then any two volume elements would be uncorrelated, which 
means that £ would vanish and the average number of pairs would simply be 
the product of the average amount of objects in each volume element. If £ 
is positive (negative), then we say those two volume elements are correlated 
(anti-correlated). 

Assuming a statistically homogeneous Universe, £ can only depend on the 
difference vector r\ 2 = r 2 — r\. Further assuming statistical isotropy allows 
£ to only depend on the distance r between 7*1 and r 2 . Hence, we denote the 
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two-point correlation function as £(r). Solving eq. (2.27) leads to 

(nin 2 ) 



- 1 = (5m(ri)5 m (r* 2 )}, 



(2.28) 



which can be easily checked by plugging in the definition of the density 
contrast 

rii - (m) 



dr. 



(2.29) 



Thus, the correlation function is often written as the average over all possible 
positions 



i y 5 m (cc)<5 m (cc + r) d 3 :r. 



(2.30) 



Another important tool that will later prove to be invaluable is the power 
spectrum, which is in a cosmological context the square of the Fourier trans- 
form of a perturbation variable (up to some normalization constant). For 
instance, the matter power spectrum is defined as 



P m (k) = V6 m (k)6* m (k) 



V 



5 m (x)e 



-ik-x j3 



d d x 



(2.31) 



Rewriting the norm yields 



Pm(k) 



1 

V 



5 m (x)5 m (y)e- lHx - y) d 3 xd 3 



or, if we substitute r = x — y 



Pm(k) 



1 

V 



S m (x)S m (x - r)e ikr d 3 x d 3 r 



-ikr j3 



d r 



(2.32) 



(2.33) 



where we used eq. (2.30) in the last step. Hence the power spectrum is the 
Fourier transform of the two-point correlation function. 

The power spectrum is amongst the cosmologist's favorite tool to describe 
our Universe in a meaningful way. It is often applied for instance to the 
matter distribution, the anisotropies of the cosmic microwave background 
radiation or the shear field of weak lensing. 

The convergence power spectrum describing the cosmic shear of galaxy 
images can be expressed in terms of the matter power spectrum. Since the 
matter distribution is 3-dimensional, but the convergence is a function of the 
2-dimensional sky, we need Limber's theorem to relate those two. It states 



4 This is the Wicner-Khinchin Theorem. 
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that the power spectrum for a projection 



F(9 X , 9y 



w(r)f(9 x r, 9 y r, r) dr . 



(2.34) 



where w(r) is a weight function (normalized to unity), turns out to be 



P(q) 



dr 



w{rY (q 

-v 



o 



(2.35) 



where p(k) is the power spectrum of /. This theorem is directly applicable to 
the total convergence field in eq. (2.26), so that its power spectrum becomes 

"OO 



4„ 

-y foe 



dr — P^ 











W(zl p (q 



with 



W(z) = 



w(r(z)) 
r(z) 



(2.36) 
(2.37) 

(2.38) 



being the window function. All we need to know now is the power spectrum 
of i[) t ij, which is simply 

P+ M = ^k^\ 2 (2.39) 

because the Fourier transform of ip^j is —kikjip. Then we can plug in the 
Poisson equation in Fourier space, which is 

k 2 ^ = 3a 2 H(a) 2 n(a)5 m (a) (2.40) 

(in the absence of anisotropic stress, i.e. ip = 2$) to obtain 

P^Jk) = k^ = m{a)^ 2 J{l + zfP m {k) . (2.41) 

Putting it all together, and replacing q with the multipole £/n, finally yields 
the power spectrum for the convergence, 



P(£) 



9H 3 



W{z) 2 E{zfn rn {z) 



-P 

J. yy 



irr(z) 



dz. 



(2.42) 



More details can be found in A&T (ch. 4.11, 14.4) or Hu and Jain (2004). 



In fig. |2.6| we can see the matter power spectrum derived from linear 
perturbation theory as well as non-linear corrections. Fig. |2.6| shows the 
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convergence power spectrum, derived from the linear and non-linear matter 



power spectrum according to eq. (2.42). 



2.4 Fisher matrix formalism 



The last tool we will need is the powerful theory of Bayesian statistics. In 
physics we define a theory by a set of parameters encapsulated in the vector 
0. The ultimate goal is to deduce the value of by performing a series 
of experiments that yield the data x. If we were given the values of 0, 
we could calculate the probability density function f(x,0) and thus predict 
the probability P of getting a realization of the data x given the theory 0, 
which is usually denoted by P(x\6). But since we rarely are given the exact 
theory, all we can do is draw samples Xi of an unknown PDF by conducting 
experiments in order to estimate the parameters 0. In a frequentist approach, 
we would take those parameters as the true ones that maximize the so-called 
likelihood function 

CfoxJ^llfixuO), (2.43) 

i 

which is nothing but the joined PDF, with now being interpreted as a 
variable and fixed parameter. But in cosmology, we often have prior 

knowledge about the parameters from other observations or theories that we 
would like to account for. For this, we need Bayes' theorem: 

s P(0\I)P(x\0,I) 

p w*^= 41/) • < 2 ' 44 > 

In this Bayesian approach, we turn the situation around so that we are 
asking the question: What is the probability of the theory, given the observed 
data? Note that we included the background information I in every term 
even though it is hardly relevant to the derivation. It is to remind ourselves 
that we always assume some sort of prior knowledge, even when we do not 
realize it. For instance, we might implicitly assume that the Copernican 
principle holds true, or that the Universe is described by the Einstein field 
equations. 

Here the prior knowledge (or just prior) is given by P(0\I), and P(x\I) is 
the marginal probability of the data x which acts as a normalization factor, 
requiring that / C(0,Xi)d n 6 = 1. P(x\0,I) is identical to the likelihood, 
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Figure 2.6 Top panel: Comparison between a matter power spec- 
trum with (purple) and without (blue) non-linear corrections in 
arbitrary units. The fitting formulae were taken from Eisenstein 



and Hu 


(1999 


) and 


Smith et al. 


(2003) 



come important only at wave numbers k > 0.1 Mpc -1 . Note that 
k is not in units of /iMpc" 1 in this case. Bottom panel: Compari- 
son between a convergence power spectrum from eq. (2.42) based 
on the linear and non-linear matter power spectrum (no redshift 
binning) . 
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which is given by the model. This could be for instance a Gaussian PDF 



n 



i=i 



er,V27r 



exp 



{Xi-Xi{9)f 
2a? 



(2.45) 



The left hand side of eq. (2.44) is then called the posterior (Hobson et al. 
2009, p. 42). We now have all the ingredients to calculate the posterior 
probability of the parameters given the data. This quantity is important for 
computing the so-called evidence, which is defined by the integral over the 



numerator in eq. (2.44) and needed for model selection. 



Unfortunately, the likelihood is often not a simple Gaussian. Finding the 



values of 6 that maximize the PDF given by eq. (2.44) can be computation- 
ally demanding, since a (naive) algorithm for finding the extremum scales 
exponentially with the number of dimensions of the parameter space. This 
will becomes a problem when considering 30 parameters as we will later in 
this thesis. It is not unusual for Bayesian statistics to suffer from the "curse 
of dimensionality" ( Bellman 2003). One way to overcome this curse is by 



using Markov chain Monte Carlo methods. In the case of a likelihood that 
is not a multivariate Gaussian, we can simply assume that it can at least 
be approximated by one. This is not an unreasonable assumption, since the 
logarithm of the likelihood can always be Taylor expanded up to the second 
order around its maximum, denoted by 6: 



ln£(0) « ln£(0) + 



d 2 ln£(0) 



dOidOj 



{(>, oniOj o,). 



(2.46) 



0=0 



No first order terms appear since they vanish in a maximum by definition. It 
is now quite useful to define the Fisher matrix as the negative of the Hessian 
of the likelihood, 

_ _ <9 2 ln£(0) 



(2.47) 



0=0 



This definition is not very useful when we need to find the maximum of 
a multi-dimensional function and then compute the derivatives numerically 
anyway, since many expensive evaluations is what we wanted to avoid in the 
first place. However, since this thesis deals with constraints by future surveys, 
using some fiducial model which we already know the expected outcome of 
and thus the peaks of the likelihood. This is the reason why the Fisher matrix 
is a perfect tool in assessing the merit of future cosmological surveys. 

Let us consider a survey in which we measure some set of observables 
x (along with their respective standard errors cr) whose theoretical values 
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x{0) depend on our model. These observables could be the same quantity 
at different redshifts, i.e. Xi = x(zi), or different quantities all together, e.g. 
x = (^cmb, c s , flk), or a mixture of both. Assuming a Gaussian PDF we can 
write down the likelihood as 



.. , I ^ (Xi - Xi(0)) 2 

C{0) oc exp I --2^- 



(2.48) 



Using the definition in eq. (2.47) we get for the Fisher matrix 

1 d 2 x k (6) 



E 



, 4 MM 



(2.49) 



0=0 



2.4.1 Calculation rules 

We will not go into a full derivation of the rules, as they can be found 
in A&T (ch. 13.3) and are quite straight forward, but rather simply state 
them here. 



Fixing a parameter. If we want to know what the Fisher matrix would be 
given that we knew one particular parameter 9i precisely, we simply remove 
the z-th row and column of the Fisher matrix. 



Marginalizing over a parameter. If, on the other hand, we want to 
disregard a particular parameter 9i, we remove z-th row and column from 
the inverse of the Fisher matrix (the so-called correlation matrix) and invert 
again afterwards. If we are only interested in exactly one parameter 6^, then 
we cross out all other rows and columns until the correlation matrix only has 
one entry left. Thus, we arrive at the important result 

a(9 i ) 2 ^(F-%. (2.50) 

This implies that the Fisher matrix must be positive definite, as it must be 
as the Hessian in a maximum. 



Combination of Fisher matrices. Including priors in our analysis is 
extremely simple in the Fisher matrix formalism, since all we need to do is 
add the Fisher matrix: 

F' = F + F prior . (2.51) 
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This only holds if both matrices have been calculated with the same under- 
lying fiducial model, i.e. the maximum likelihood is the same. The only 
difficulty lies in ensuring that the parameters of the model belonging to each 
matrix are identical and line up properly. If one matrix covers additional 
parameters, then the other matrix must be extended with rows and columns 
of zeros accordingly. 

Parameter transformation. Often a particular parameterization of a 
model is not unique and there exists a transformation of parameters 



which might happen when combining Fisher matrices from different sources. 
Then the Fisher matrix transforms like a tensor, i.e. 



& = 0'{0) , 



(2.52) 



F' = 



J T FJ, 



(2.53) 



where J is the Jacobian matrix 




d6i 



(2.54) 



which does not necessarily need to be a square matrix. If it is not, however, 
note that the new Fisher matrix will be degenerate and using it only makes 
sense when combining it with at least one matrix from a different source. 



Chapter 3 



Constraints on cosmological 
parameters by future dark 
energy surveys 

With the knowledge of the observational specifications of future weak lensing 
surveys, we can use the Fisher matrix formalism to forecast the errors that 
said surveys will place on cosmological parameters. In particular, we want to 
expand the list of those parameters with values of the Hubble parameter and 
growth function at certain redshifts. An invaluable resource for this chapter 
is the book by A&T (ch. 14) 



3.1 Power spectrum tomography 



It has been shown that dividing up the distribution of lensed galaxies into 
redshift bins and measuring the convergence power spectrum in each of these 
bins as well as their cross-correlation can increase the amount of information 



extracted from weak lensing surveys (Hu, 1999; Huterer, 2002). This means 



on one hand that we need to have additional redshift information on the 
lensed galaxies, but on the other that we are possibly rewarded with gaining 
knowledge about the evolution of dark energy parameters. Note that the 
purpose of the redshift bins in this thesis is two-fold, and power spectrum 
tomography is only one of them. The other is the linear interpolation of 
H(z) and G(z), where the centers of the redshift bins act as supporting 
points, making our analysis independent of assumptions about the growth 
function by any particular model. We choose to divide the redshift space 
into M bins such that each redshift bin contains roughly the same amount 
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of galaxies according to the galaxy density function n(z), i.e. 

-y poo 



n(z) dz = — -. / n(z) dz 
^ Jo 



(3.1) 



for each i. We infer the values of the Zi via a series of successive numerical 
integrations with zo = and = 3 (see fig. reffig:2-zbins). A common 
parametrization of n(z) is 



n(z) oc z a e (z/2p) 



(3.2) 



with a = 2 and (3 = | . Here z p is related to the median redshift z m = lA12z p 
(Amara and Refregier, 2007). We need a galaxy density function for each 



0.15 



0.10 



N 



0.05 



0.00^- 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Figure 3.1 The area under the curve of the galaxy density func- 
tion is the same in each redshift bin for N = 5. 

bin, and the naive choice would be 

~ / s / n(z) Zi_i < z < Zi, 

n ^ = { else. (3 ' 3) 

But to account for redshift measurement uncertainties, we will convolve n(z) 
with the probability distribution of the measured redshift z p h given the real 
value z: 

ni (z)= / dz'n i (z')P i (z ph = z'\z). (3.4) 
Jo 

Note that the integral vanishes if z' lies outside the ith bin, so the limits 
can be adjusted to the according finite values. This convolution reflects the 
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fact that we cannot assign a galaxy to a particular bin by measuring the 
photometric redshift with finite precision. The probability distribution is 
here modeled by Gaussian, i.e. 



Pi{z ph \z) 



1 



27TCT, 



exp 



Q P h - zY 

2a? 



(3.5) 



with <jj = A z (l + (zi_i + Zi)/2) (see Ma et al. (2006) for details). Finally, we 
normalize each function to unity, such that 



nAz) dz 



(3.6) 



The resulting functions can be visualized as in fig. 3.2 




Figure 3.2 Normalized galaxy density functions for each bin with 
N = 5, convolved with a Gaussian to account for photometric 
redshift measurement errors. 



3.2 The matter power spectrum 



Future weak lensing surveys will supply us with the convergence power spec- 
trum in which the cosmological information we are seeking is imprinted. 
However, for now we will have to rely on simulated data. As we will see in 



section 3.3, the convergence power spectrum depends on the matter power 



spectrum which can be simulated by a fitting formula derived by Eisenstein 



and Hu (1999). Non-linear corrections need to be accounted for, since we are 
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assuming an angular galaxy density of 35arcmin 2 , which corresponds to a 
multipole up to a order of magnitude of 



lg I V 35 arcmin 



(3.7) 



As we saw in fig. 2.6 non-linear corrections are required for i > 10 3 , and for 



this we need to use the results by Smith et al. (2003). Here we will give a brief 



overview of the essential results found in the two papers cited above, using 
our notation and some simplifications (no massive neutrinos, flat cosmology). 



3.2.1 Fitting formulas 



It should be stressed that the growth function D(z) in Eisenstein's paper 
differs from ours by a factor of a = 1/(1 + z), so D(z) = aG(z). A precise 
definition of the growth function will be given in section. |3.6.2 They also de- 
fine 02.7 = Tcmb/2.7K. Only in this section will we differentiate between the 
linear and the non-linear power spectrum. Later on, the expression "power 
spectrum" or P m (k) will always imply that non-linear corrections have been 
applied. 

Using the definition of the growth function, the matter power spectrum 



in the linear regime can be written as (cf. eq. (3.70)) 



1 + zY 



(3.8) 



where n s is called the scalar spectral index. The transfer function T(k) can 
now be fitted by a series of functions as follows. 



T(k) 



LB(k) 
L + Cg c 2 ff ' 



(3.9) 



with (we assume that there are no massive neutrinos, in which case B(k) 
equals unity) 



L = ln(e + lMPcy/c^qes) 
67= 1.44+ 325 



1 + 60.5^ n 



(3.10) 
(3.11) 
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where 



q c s 



kQ 2 



2.7 



r cfr M pc - ] 

r cfr = u m | yja v + 



1 - Ja^ 



1 + (0.34ibs)' 



and 



a 



v = — 5 2(PC / Pm) (1 - 0.553u; 6 + 0.126^) (1 + y d f^ 

1 

1 + 



x <M + 



'm 

Pc ~ Pcb 



(3 - A Pc )(7 - Ap m ) 



Here, we used the abbreviations 



p x = -(5- Vl + 24w x ) 



Vd = 



l + z, 



cq 



with 



z ccl = 2.50 x lO 4 w m 2 .7 



z d = 1291 



,0.251 



(3.12) 

(3.13) 
(3.14) 



(3.15) 



1 + 0.659c^ 828 

61 = 0.313w-°- 419 (l + 0.607^ 674 ) 

6 2 = 0.238w^ 223 . 

By definition, the power spectrum needs to be normalized such that 

r0 ° dA; k 3 



4= / 



k 2tt 2 



P L (fc)|H/ 8 (A;)| 2 , 



(3.16) 
(3.17) 

(3.18) 

(3.19) 

(3.20) 
(3.21) 

(3.22) 



where Wn(k) is the Fourier transform of the real-space window function, in 



40 



Chapter 3. Constraints by future surveys 



this case a spherical top hat of radius R (in Mpc), i.e. 

W «M K { I ottaJL < 3 - 23 > 

3 

W R (k) = — — ^(sin kR- kR cos kR) . (3.24) 
3.2.2 Non-linear corrections 

When applying non-linear corrections, we use the dimensionless power spec- 
trum, which is defined by 

Al(k) = —P h (k) , (3.25) 

and similarly for other indices. It turns out that the non-linear power spec- 
trum can be decomposed into a quasi-linear term and contributions from 
self-correlations 

A* L = A 2 Q (k) + A'(fc) , (3.26) 

which are given by 

A «w = A '^» + exp( " s/4 " yV8) (3 - 27) 

with y = k/k a and 

A 2 ' _ ^ny_ (o nq) 

H 1+6ny / 2 (n m) + (Cn/ 3 (fim)y) 3- 7n - ^ • ^ 

In these equations, k c is defined via 

a(R G ) 2 = A 2 L (k)&cp(-k 2 R G ) 2 dlnk (3.30) 

J oo 

by 

v(K 1 ) = l, (3.31) 

while the effective index is 

din a(R) 2 



n = — 3 — 



dlni? 



(3.32) 



(7=1 
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and the spectral curvature is 

d 2 In a {Rf 



C 



dlni? 

The best fit yielded the following values for the coefficients: 



(3.33) 







= 1.4861 + 1.8369n + 1.6762n 2 + 0.7940n 3 + 0.167n 4 


- 0.6206C 


k 




= 0.9463 + 0.9466n + 0.3084n 2 - 0.9400C 


(3.34) 


k 




= -0.2807 + 0.6669n + 0.3214n 2 - 0.0793C 


(3.35) 




o; n 


= 1.3884 + 0.3700n - 0.1452n 2 


(3.36) 




fin 


= 0.8291 + 0.9854n + 0.3401™ 2 


(3.37) 




In 


= 0.8649 + 0.2989n + 0.1631C 


(3.38) 


Ig 




= -3.5442 + 0.908n 


(3.39) 


k 




= 0.9589 + 1.2857n. 


(3.40) 


Also, 


the 


functions ft in a flat universe are given by 










(3.41) 






f2{n m ) = fi- - 0585 


(3.42) 






/ 3 (fi m ) = fir 43 - 


(3.43) 



3.2.3 Parameterized post-Newtonian formalism 

We shall allow another degree of freedom in our equations that stems from 
scalar theories in more than four dimensions. The Gauss-Bonnet theorem in 
differential geometry connects the Euler characteristic of a two-dimensional 
surface with the integral over its curvature. In general relativity, it gives 
rise to a term with unique properties: It is the most general term that, when 
added to the Einstein-Hilbert action in more than four dimensions, leaves the 
field equation a second order differential equation. In four dimensions, the 
equation of motion does not change at all under this generalization, unless 
we are working in the context of a scalar-tensor theory, in which the Gauss- 
Bonnet term couples to the scalar part and modifies the equation of motion. 
For our purposes, this will effectively make Newton's gravitational constant 



G a variable, denoted by G*{rj) (for details, see Amendola et al. (2006) and 
references therein). This is called the parameterized post-Newtonian (PPN) 
formalism. 
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By defining 

Q = ^r, (3-44) 
one can write the Poisson equation in Fourier space as 

k 2 <f) = -A7cGa 2 Q Pm A m , (3.45) 

where A m accounts for comoving density perturbations of matter. If we 
admit anisotropic stress, then the two scalar gravitational potentials do not 
satisfy ^ = — $, and we parameterize this by 

= -(1 +r](k,a))$. (3.46) 

With weak lensing, only the combination 

£ ee Q(l + 77/2) (3.47) 

appears in the convergence power spectrum in a very simple way, where we 
can just replace the matter power spectrum P m (k) by HP m (k). Obviously, 
in the standard ACDM model, we have £ = 1, so we allow it to vary in time 
and parameterize it as 

S(a) = l + S a, (3.48) 

such that initially it looks like the standard model and then gradually di- 
verges from there (Amendola et al, 2008). Finally, we add E to our list of 
cosmological parameters with a fiducial value of 0. 



3.3 The Fisher matrix for the convergence 
power spectrum 

Building on the methods outlined in chapter 2, we can now compute the 
weak lensing Fisher matrix. A full derivation of its expression can be found 



in A&T, ch. 14.4 or Hu (1999), with the final result being 



^ (21+ l)Al ^ dP,(£) . dPU?) r -i ( ~ AQ s 

V - /sky 2^ 7) 2^ Q0 q C jk Qf L mi ■ (3.49) 



ijkm 



For this we need a multitude of cosmological functions that will be defined 
in this section from the bottom up. First of all, note that the sum should 
go over all multipoles £ inside an interval £ m in--^max that is determined by 
the fractional survey size / s k y and the angular galaxy density n$. However, 
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this poses a computational challenge. Hence, we rather sum over bands with 
width A£ while keeping in mind that one multipole £ contains (2£ + 1) modes 



(Hu and Jain, 2004). The multipole intervals will be logarithmically spaced 



and the choice lg A£ shall be justified in section 4.2.2 
Cjk is the covariance matrix, given by 

C jk = P jk + Sjk^nJ 1 , 



(3.50) 



where 7; n t is the shot noise due to the intrinsic ellipticity and rij is the number 
of galaxies per steradians belonging to the j-th bin: 



Tli 



3600 



180 



7T 



n$ / rij(z)dz, 



(3.51) 



where n$ is the galaxy density per arcmin 2 . The convergence spectrum Py 
depends on the matter power spectrum P m and is given by (now with the 
PPN parameter) 



P m _ 9 f 00 W^W^H^z^zf 



dz 



(3.52) 



which can be slightly simplified, by using Q m (z) = Q m x (1 + z) 3 /E(z) 2 , to 

WiWWjWil + z) 2 



Pi, 



A lL m 



E(z) 



~^{z)P„ 



irr(z) 



dz. (3.53) 



Here, H(z) denotes the Hubble parameter, which is given by the first Fried- 

(3.54) 



mann equation (eq. (2.1)) as 



H(z) = # v^Ul + zf + (1 - Q m ) exp(f w (z)) , 
and E(z) = H(z)/H is its dimensionless equivalent. The function 

l + w(z) 



dz'. 



(3.55) 



describes how the dark energy density scales with the scale factor and is the 



solution of the continuity equation (2.5). The equation-of-state ratio w(z 



p/p is from now on assumed to take the common CPL parameterization 1 



Relevant literature often uses w a instead of w\. 
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QChevallier and Polarski] [200lj pnd"er] [2003]) 



w(z) = W + W\ 



1 + z 



(3.56) 



in which case the above function can be calculated analytically and takes the 
form 

W\Z 



f w (z) = 3[(w + w 1 + l)ln(z + l) 



z + l 



(3.57) 



One of the most crucial functions being used in this calculation is the window 




function 



WAz) 



dz' 



1-4^-1 n,(r(/)), 



H(z') \~ r(z') 
where r(z) is the comoving distance at redshift z, i.e 

r(z) 



H(z' 



(3.58) 



(3.59) 



Lastly, we need the matter power spectrum P m to infer the convergence power 



spectrum. For this we use the fitting formulae covered in section 3.2 



3.4 Survey parameters 



45 



3.4 Survey parameters 

The quality of any weak lensing survey is determined by a set of parameters 
that includes the fraction of the sky covered 2 , the median redshift, and so 



on. In table 3T all of the parameters that are used in our simulation are 
shown. Weak lensing surveys should aim at wide and deep field coverage. 
Since these projects are all still in development, the exact numbers that will 
be used in the end may vary, and some are still undecided. 



Parameter 


Value 


Comment 


/sky 


0.5 


Fraction of sky surveyed 




0.9 


Median redshift 


A, 


0.05 


Relative photometric redshift error 




35 


Galaxy density per arcmin 2 


Tint 


0.22 


Intrinsic shear 



Table 3.1 Shown here are typical values of survey parameters 
determining the quality of the experiment, which we will use in 



our calculation (cf. Huterer et al. , 2006 ) 



3.4.1 Ground based surveys: LSST 

The Large Synoptic Survey Telescope, 3 funded partially privately and by 
the U.S. Department of Energy, is a wide field survey telescope working in 
the visible band, located on the Cerro Pachon in Chile. It is currently in the 
design and development stage and is scheduled to be fully operational by 2020 
with a lifetime of at least ten years. It will feature a 3.2 Gigapixel camera 
and a large aperture with an 8.4m (6.5m effective) primary mirror and cover 
the entire extra galactic sky of 20 000 deg 2 with a depth up to an apparent 
magnitude of r ~ 27.7 for point sources. Besides weak lensing measurements, 
the mission will observe super novae as well as map the Milky Way and small 
objects in the solar system (such as asteroids). This ultra wide and deep 
survey will provide a good data source for weak lensing observations. 



2 Since the Milky Way is blocking roughly half of our sky for deep surveys, even the 
most ambitious missions will have / s k y = 0.5. 
http://www.lsst. org/ 
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3.4.2 Space based surveys: Euclid, WFIRST 



Euclid 4 is a proposed mission from ESA for a satellite orbiting the second 
Lagrangian point as part of their cosmic vision program. It is still in the 
definition phase and currently competing with PLATO and Solar Orbiter to 
be one the two missions that will be approved by ESA in mid 2011. It is 
foreseen to be launched in 2017 with a nominal mission lifetime of five years. 
This survey might also cover up to 20 000 deg 2 . The telescope will be a 1.2m 
Korsch operating within the visible and infrared wavelengths with the ability 
of detecting galaxies in the redshift range 0.5 < z < 2 and therefore provide 
an excellent probe for dark energy. 

The Wide Field Infrared Survey Telescope (WFIRST 5 , formerly known 
as JDEM-Omega) is another survey investigating dark energy, funded by 
NASA and the U.S. Department of Energy and proposed by the Joint Dark 
Energy Mission. In addition to weak lensing and baryon acoustic oscillations, 
WFIRST will also probe supernovae with high precision. It is still at a very 
early stage and has been postponed to launch at an unknown date, possibly 
not before 2025, due to budgeting and scheduling issues on a seperate project, 



the James Webb Space Telescope (Overbye, Jan 4, 2011). 



Survey f t 



sky 



A, 



n$/ arcmin Reference 



LSST 0.37 0.02 0.7 

Euclid 0.5 0.05 0.8 
WFIRST 0.25 0.04 ? 



30 

30-40 
> 30 



Huterer et al. (2006); 



Ivezic et al. (2008) 



Laureijs| (|2009) 
Gehrelslteuf 



Table 3.2 Overview of current and planned surveys. 



3.5 Fiducial model 

In order to apply the Fisher matrix formalism, we need a fiducial cosmological 
model which represents the vector in parameter space where the likelihood 
has its maximum. We use the standard model of cosmology, a flat universe 
with cold dark matter and a cosmological constant (called ACDM, which 
makes six parameters). We additionally consider some extra parameters 
which are allowed to vary in a non-standard way, in particular 

= (u m , u b} r, n s , Q m , w , w h 7, S , er 8 ) , (3.60) 

A — 1 

http : //sci . esa. int/science-e/www/object/ index . cfm?f object id=42266 
5 http : //wf irst . gsf c . nasa. gov 
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Parameter 


Fiducial value 


Description 


UJ b 
Q m 


0.1352(36) 

0.02255(54) 

0.275(11) 


Reduced matter density today 
Reduced baryon density today 
Matter density today 


T 


0.088(14) 


Reionization optical depth 


n s 


0.968(12) 


Scalar spectral index 




0.816(24) 


Fluctuation amplitude at 8Mpc//i 


wl 


-1 


Current equation-of-state ratio 


w\ 





Higher order equation-of-state ratio 


7* 


6/11 


Growth index 







Parameterized Post-Newtonian parameter 



Table 3.3 Our ten-parameter fiducial model. The second col- 
umn represents the WMAP7+BAO+-£fo Mean taken from Ko- 



matsu et al. (2010). Quantities with a dagger are fixed in the 



standard flat ACDM model and thus not included here by the 
WMAP measurement. Uncertainties in the last two digits are 
given in parenthesis. They will be needed in Section 4.5 



where u}„ 



is the reduced fractional matter density, and Ub is defined 



analogous for baryons. Their numerical values according to the WMAP 7- 

The growth index 



year data (Komatsu et al. , 2010) can be seen in table 



7 is taken to be 6/11 ~ 0.55 in the standard model (Wang and Steinhardt 



3.3 



1998). 



It should be pointed out that the reionization optical depth r is merely 
a remnant, and none of the functions in the code in later section depend 
on r; however, it was not removed, however, due to the possibility of future 
modifications to the code. The inclusion of functions that depend on r will 
thus be easier, as long as we remember to cross out the according row and 
column in the resulting Fisher matrix. 



3.6 Model-independent parameterization 
3.6.1 Expansion rate 

In order to see how observations can constrain the history of the expansion 
rate without assuming a particular parameterization of it, we take the values 



of the Hubble parameter at a series of redshifts determined in section 3.1 



using our fiducial model and build a linearly interpolated function going 
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through these supporting points. For this, we define 

fH = ]nH(zi), (3.61) 

where Zi is at the center of the z-th redshift bin. The use of the logarithm 
is purely for convenience. The Hubble parameter itself from now on is a 
linearly interpolated function through these points. It is evident that the 
dependency of the Hubble parameter on the other cosmological quantities 
changes from H(z; fl m , w , W\) to H(z; hi, ... , hj\f). The consequence of this 
is that all cosmological functions that depend on H(z) are now treated as 
functions of the his, resulting in M additional parameters being added to the 
fiducial model. 



3.6.2 Growth function 

The growth function is defined by 

G{a) 



$(a T ) 



(3.62) 



where $ is the scalar gravitational potential from (2.8) and ax stands for the 
transfer epoch, which typically has a value of 0.03 (see for instance A&T or 



Dodelson (2003)). Be advised that in relevant literature the growth function 



is often defined such that G(a) needs to be replaced with aG(a) in eq. (3.62). 
Also, it is sometimes referred to as D or D\. The growth function can be 
identified with the matter density contrast over the scale factor, i.e. 



G(a) 



which is given by 



a = 



a5 m (0) 

pip 



(p) 



- 1 



(3.63) 



(3.64) 



A practical quantity in cosmology is the so called growth index 7, defined 



by 



where / is the growth rate 



f(a) = Vt m (a)' 



f(a) 



dln<L 



din a 



(3.65) 



(3.66) 



Note that the power law in eq. (3.65) is an empirically found fit for the 
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growth rate (Wang and Steinhardt, 1998), which was shown to be accurate 



to within 0.2% of the exact solution (Linder, 2005). 
After reminding ourselves that 



din a 



da 



[l + z)d 



1 



1 + z 



dz 
1 + z' 



(3.67) 



we can eliminate 5 m (a) and / from eqs. ( |3.62| - |3T65| ) to find 

G(z) = Go exp 
where the matter density is 



q»(^) 7 - 1 

l + z' 



dz' 



(3.68) 



Pm{z) 



Pm(z) +PA 



Q m x (1 + z) 
E(z) 



(3.69) 



The only point where the growth function enters the calculation is via 



the matter power spectrum, which is given by (see A&T, p. 75; cf. eq. (3.8)) 



Pm\k, z) 



2n 2 5 2 H ( k 



,2 G(zf 



H J T ^lffW H ° 



Just as in section 3.6.1 we can now finally define 

gi = InG(zi) 



(3.70) 



(3.71) 



and replace the original growth function with a linear interpolation through 
the 2iS. The number of parameters in our model thus increases to 10 + 2j\f. 
Naturally, every function that depends on the matter power spectrum now 
also becomes a function of the ^s. 



3.7 The figure of merit 



In order to have a quantifiable gauge of the information content of the final 
result, we define the figure of merit (FOM) for the error bars on the h^s and 
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<7iS, respectively, as 



FOM(#) = v( h i)~ 2 > (3-72) 



i=l 

N 

FOM(C7)^^rr(^)- 2 . (3.73) 
i=i 

It is obvious that the smaller the error bars are, the larger the FOM becomes. 
The idea of using this definition is that the FOM should not change too much 
if only one error bar is significantly larger than all the other ones. 

While the particular value of the FOM has no deeper meaning, it provides 
a convenient tool to easily compare the errors resulting from different sets of 
parameters. Note that this definition does not match the popular definition 



by the Dark Energy Task Force, which suits a different purpose (Albrecht 
eToO [20061). 



3.8 An implementation in Mathematica 

To do the actual computation, we use Mathematica 8.0.1.0, which supplies 
us with all necessary numerical procedures. 

In an attempt to seperate all the parameters that determine the final 
outcome of the computation from the rest of the code, they are all defined 
at the beginning of the Mathematica notebook. We further calculate a 
checksum, a digital fingerprint of those parameters via an injective 6 map 
from the set of all parameters onto the set of integers, comparable to a hash 
function. This will be useful when we store results in a file when changing 
the parameters without having to worry about overwriting previous results. 
We can roughly differentiate between physical parameters (such as matter 
density, baryon density, etc.) and technical parameters (the limit of a certain 
integral, step size in an interpolating function, etc.). For an overview, we refer 



to tables 3.4, 3.5, and 3.6 



In the first section of actual computation, we infer the boundaries of the 
redshift bins depending only on Af, A z , and the galaxy distribution func- 
tion n(z). The definition of the necessary cosmological functions is straight 
forward. They often involve the computation of time intensive numerical 
integrals, especially the highly used window function. These functions are 
interpolated by cubic splines. All functions are sufficiently "nice" (not more 



3 At least for our purposes. 
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Functions 



logprint [string] 
loadf ile [filename] 
ngal [ [i] ] [z] 
hub [z, dbins] 
dist [z,dbins] 
omgen [z , om , wO , wl] 
growthf astgen [z , gbins] 
wind [i , z , om , dbins] 
pk[k,z, . . .] 
pknl [k , z , . . . ] 
clnngen[i, j ,ell, . . .] 
dcdp [pi , i , j , ell , ref ] 



Prints messages and saves them to a log file 
Uses pre-calculated expressions if applicaple 

rii(z) 

H(z;h 1 ,...,h Af ) 

r(z) 

Vt m (z) 

G(z;gi,...,gx) 
Wi(z) 

P m (k) (Linear) 
P m (k) (Non-linear) 

PM 



Table 3.4 Functions as used in the Mathematica notebook. 



Physical parameters 



omh2 


L0 m 


obh2 




tau 


T 


ns 


n s 


om 




wO 


w 


wl 


Wi 


gamma 


7 


gPPn 


So 


s8 




dbins [[i]] 


hi 


gbins [[i] ] 


9i 



Table 3.5 Variables referring to physical parameters as used for 
arguments in functions in the Mathematica notebook. Their 
fiducial value is stored in a variable that has the string ref at- 
tached. 



than one extremum) so that we do not have to retreat to more conservative 
linear splines. 

In order to save processing time and memory (especially in sequenced ex- 
ecutions of the notebook) we employ a little trick when computing expensive 
functions. To avoid loss of data after unexpected computer crashes, it makes 
sense to store not only the final result, i.e. the Fisher matrix, but also inter- 
mediate results on the hard drive. For this, the function loadf ile evaluates 
the expression given in the second argument only if the file given in the first 
argument does not exist, and then stores the result in said file and returns it. 
If the file does exist, the expression in it will be loaded and returned. This 
is applied to all time sensitive functions. All files are saved in dir, which 
contains the checksum of all parameters, except for the window functions 
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Technical parameters 



zlimit 


3 


Redshift limit after which we assume n(z) ~ 


nbins 


5 


Number of redshift bins J\f 


eps 


0.025 


Step size for numerical differentiation 


zsigma 


4 


Integration length in units of standard deviations 


zsmax 


4 


Maximum redshift for interpolated functions 


f itstepz 


0.2 


Step size for interpolated functions 


lellmax 


4 


Decadic logarithm of upper limit of the multipole 


lellmin 


2 


Decadic logarithm of lower limit of the multipole 


lstep 


0.2 


Decadic logarithmic step size of the multipole 


piv 




List of parameters that are not fixed 


dir 




Directory where data is stored 



Table 3.6 Parameters determining the accuracy of the result and 
their default value. 



Wi(z). Because they only depend on fl m , J\f, the h^s and the number of the 
redshift bin, they are stored in a sibling directory to save computing time 
when calculating the Fisher matrix after only, say, A lg I has changed. The 
computation of the interpolating functions of the window functions has also 
been parallelized by using the M ATHEMATICA command ParallelTable. 

Very often we have to interpolate a function depending on several pa- 
rameters, such as the window function Wi(z; fl m , hi, . . . hjj), which generally 
needs to be evaluated for the derivatives at all < z < 3 for Q m , fi m (l ± e) 
and so on. We define the function such that it is computed only when needed, 
which can be achieved in Mathematica by writing schematically 

fAux[a_] := fAux[a] = 

Interpolation [Table [{x,f [x,a] }, {x,xO,xl ,dx}] ] ; 
fFast[x,a] := fAux[a] [x] ; 

The fFast[x,a] is a fast, interpolated version of f [x,a], where x is the 
variable and a a parameter. The extra = sign in the first line makes Math- 
ematica look up the cached value if there is one, preventing it from making 
an expensive computation more than once. 



What follows is a relatively straight forward implementation of sections 3.2 



and 3J3 It is important to note, however, that all the functions that ulti- 
mately depend on the growth function, i.e. all functions that depend on the 
linear matter power spectrum, have been overloaded such that they might 
take gbins as an argument, in which case the power spectrum uses the in- 
terpolated growth function. The reason for this will become apparent in 
chapter 4. 
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Most functions have been carefully defined to avoid warnings by Mathe- 
matica concerning numerical issues like non-converging integrals. Vanishing 
integrals are hard to detect with numerical methods, since a sufficiently sharp, 
high peak would require an arbitrary large number of sampling points. Many 
functions have been tweaked (if possible) to avoid these warnings. Very often, 
Mathematica also tries to simplify an expression before plugging in actual 
values, which causes a lot of warnings and can slow down the computation 
significantly. To prevent this, we declared certain arguments as "numerical" 
so that Mathematica plugs in these values immediately. For instance the 
comoving distance is defined as 

dist [z_?NumericQ , dbins_] := NIntegrate [1/hub [zx, dbins] , 
{zx, 0, z}] ; 

This way, MATHEMATICA will not try to numerically evaluate the integral 
in a subsequent definition of a function (which would be futile) before a 
specific value for the redshift is given. We found that this significantly reduces 
numerical instabilities and improves the overall results. 



Chapter 4 
Results 



Before we cover the analysis of the results, it needs to be pointed out that 
two Fisher matrices have been calculated. The first Fisher matrix was com- 



puted while using the standard expression in eq. (3.68) for G(z) and linearly 
interpolating H(z), and the second Fisher matrix was computed by linearly 
interpolating both H(z) and G(z). We will refer to this as the H-case and the 
G-case respectively. This introduces a certain degeneracy, since the variables 
hi occur in both the H-ca.se and the G-case and generally take on different 
values. Hence we will denote those variables in the if-case as hi to avoid 
confusion, although they will be barely used. 

Unless mentioned otherwise, the Fisher matrix has been calculated using 
the parameters shown in tab. 3.6[ where we cross out the rows and columns 
corresponding to (r, w , wi, 7, £ ) in the standard case. The order of our 



parameters was given in eq. (3.60); here is a reminder: 



= (u m , u b , t, n s , Q m , w Q , wt, 7, S , <x 8 ) . 



(4.1) 



From the resulting Fisher matrix, we can easily extract the errors on each 
cosmological parameter. A typical example of plotting the error bars on the 



Hubble parameter and growth function can be seen in fig. 4.1 



4.1 Consistency checks 

Since we naturally cannot compare our results with observations within the 
next decade, we have limited means to make sure that they are meaningful 
and error free. One way to check is to see whether the matrix satisfies basic 
properties, another is to compute intermediate, partial results and compare 
them with the numerical results. A first approach would be to visually com- 
pare the derivative of the convergence power spectrum in two particular bins 
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-0.5 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 

z 

Figure 4.1 A typical plot of the error bars on the his (top) and 
giS (bottom) for N = 5 with a figure of merit of 7634 and 2865 
respectively. The bins are chosen such that they all contain the 
same number of galaxies. 



with respect to a particular cosmological parameter, e.g. 

dPM 
dgi 



(4.2) 



with and without interpolated functions. However, this turned out to be 
unpractical due to the long runtime of over three days even with the use 
of four parallel kernels and reducing the sampling points in the numerical 
integration down to 50. 

But we can still easily do another check. As the negative of the Hessian 
of a function at its maximum, the Fisher matrix needs to be positive defi- 
nite. Looking at the eigenvalues of the Fisher matrix, we find in the if -case 
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(descending order) 

X H = (1.03 x 10 8 , 4.87 x 10 6 , 1.89 x 10 5 , 7.29 x 10 4 , 3.34 x 10 4 , 
9.11 x 10 3 , 6.71 x 10 3 , 6.57 x 10 2 , 3.2 x 10 2 , 1.06 x 10 2 , 

4.76 x 10 1 , 1.67 x 10 1 , 3.53, 1.04 x 1(T 2 , 1.4 x 1(T 15 ) . (4.3) 

Note that the last eigenvalue is thirteen orders of magnitude smaller than the 
next largest eigenvalue. The corresponding normalized eigenvector is, while 
neglecting numbers that are smaller than 10 -13 , 

e H ,is = (0, 0, -1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) . (4.4) 

This reveals that it points only in the r-direction in parameter space, as 



we can see in eq. (4.1). This degeneracy is expected, since the convergence 
power spectrum does not depend on r. All other eigenvalues are positive, so 
this check is passed. 

In the C7-case, the eigenvalues of the Fisher matrix are 

A G = (1.24 x 10 8 , 5.21 x 10 6 , 6.63 x 10 5 , 1.86 x 10 5 , 8.85 x 10 4 , 
4.08 x 10 4 ,9.26 x 10 3 ,7.65 x 10 3 ,2.96 x 10 3 ,1. x 10 3 , 

5.16 x 10 2 , 1.4 x 10 2 , 30.3, 7.36, 3.56, 0.83, 
- 1.92 x 1(T 14 , 5.33 x 1(T 15 , -5.31 x 1(T 15 , 1.07 x 1(T 46 ) . (4.5) 

There are two negative eigenvalues among the last four of them, but as 
they are at least thirteen orders of magnitude smaller than the next largest 
eigenvalue, they are consistent with zero. Again, analyzing the according 
eigenvectors while neglecting values below 10~ 13 yields 

e G>17 = (0, 0, 0, 0, 0, 0.085, -0.23, 0.97, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) (4.6) 

e Gjl8 = (0, 0, 0, 0, 0, 0.98, 0.17, -0.047, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) (4.7) 

e G19 = (0, 0, 0, 0, 0, -0.15, 0.96, 0.24, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) (4.8) 

e G ,2o = (0, 0, 1.0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) . (4.9) 



Comparing again with eq. (4.1) shows that these eigenvectors lie in the pa- 
rameter subspace spanned by r, w , W\ and 7, as we would expect since 
now the growth function depends on the g^s instead of the equation-of-state 
parameters and the growth index 7. At this point, the convergence power 
spectrum does not depend on any quantity that in turn depends on those 
parameters, so the Fisher matrix will be degenerate if we do not cross out 
the according rows and columns. Just like in the if-case, all remaining eigen- 
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values are positive. 



4.2 Improving the figure of merit 

In this section we want to find out how the FOM behaves when we modify 
various parameters. This provides us with a way of gauging the quality of the 
tomography process, seeing how much information we were able to extract 
out of the weak lensing signal. 



4.2.1 Impact of the number of redshift bins 

Running the MATHEMATICA notebook for various values of M lets us com- 
pare the FOM for different binning. To get optimal results, we let the number 
of redshift bins M run between 2 and 20, using the values 2, 3, ... , 10, 15 
and 20. In fig. 4J2 we can clearly see a peak at M = 3 in the H-case with a 
decreasing trend for higher numbers of redshift bins. 



12000 : 
10000 : 
8000- 
6000- 
4000 : 
2000- 

oU 



2 3 4 5 6 7 
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9 10 15 20 



Figure 4.2 Comparing the FOM for various values of M in the 
-ff-case. 



Doing the same in the G-case yields a slightly different result, shown in 
fig. O) Here the FOM for H peaks at N = 6 while the FOM for G peaks 



at Af = 3. 

When examining the galaxy distribution function for each bin (see fig. 4.3 ) 
in the case of 15 redshift bins, we can see that with our assumed photometric 
redshift error they are "smeared out" to the point where it becomes impos- 
sible to confidently assign a redshift bin to a given galaxy. 
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Figure 4.3 The galaxy distribution function in 15 redshift bins, 
which are overlapping heavily. 

4.2.2 On £ max and A\g£ 



It is not obvious how far the sum in eq 
size for £ should be. Estimating a realistic value for the upper limit £ 



(3.49) should go or what the step 

max has 

been done in eq. (3.7), but it is interesting to see how the FOM behaves for 
larger values. We therefore try several values for £ max and see whether the 



FOM approaches some constant or not. The result can be seen in fig. 4.5 



We observe that the FOM only changes slightly for values above 4, as we 
would expect due to the finite angular galaxy density. 

Doing the same for some values for A\g£ at a fixed £ max = 10 4 yields the 



numbers used in fig. |4.6| We see that essentially all the information is being 
extracted at a logarithmic multipole step size of 0.2. This figure roughly 



matches the result of 0.3 found in Bernstein (2009), where it is argued that 



the broad window function "smoothes away fine structure" in the convergence 
power spectrum, such that at some point further increasing of the number of 
multipole bins will no longer have beneficial effects. On the other hand, it is 
interesting to note that the FOM increases almost by an order of magnitude 
by going from A \g£ = 0.8 to 0.4. 



4.3 Fixing various parameters 

In order to determine the parameter which the Fisher matrix is most sen- 
sitive, we will fix one additional parameter after another and analyze the 
behavior of the FOM. Remember, fixing a parameter within the Fisher ma- 
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Figure 4.4 In the G-case, the FOM in dependence of M takes 
a slightly different shape. Blue represents FOM(H), green is 
FOM(G) and yellow is the sum of both. 
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Figure 4.5 When increasing 



x , the upper limit of the sum 
in eq. (3.49), the FOM changes, but increases by less than 15% 
beyond £ m ax = 10 4 . Shown here is the -ff-case, the G-case looks 
similar. 
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Figure 4.6 The figure of merit for various values of Alg^ in the 
H-c&se. Going beyond 0.2 does not change much. The situation 
is the same in the G-case and thus not shown here. 
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trix formalism is as easy as crossing out the according row and column of the 



Fisher matrix. As we can see in fig. 4.7 knowing the scalar spectral index 
n s in the H well as possible pays off the most, as the FOM increases 

roughly by a factor of three if we know it exactly. 



On the other hand, in the G-case, depicted in fig. 4.8, it is most beneficial 
to know Q m (factor of six) to get the highest increase in FOM(if), but 
knowing <jg exactly gives another factor of two for FOM(G). 



7000C 
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Figure 4.7 Various parameters are being fixed successively {H- 
case). Each bar is labeled with the parameter that has been fixed 
on top of the previous one, such that all of them are fixed on the 
right side of the diagram. Fixing the scalar spectral index is most 
beneficial, as the FOM triples in the best case. 



4.4 Uncertainties 

Finally, we can present the uncertainties on the values of H and G at different 
redshifts. We use the parameters Af = 5, £ max = 10 4 and A\g£ = 0.2 and find 
the results shown in tables |4~T1 and l4~2l These are derived from the full Fisher 



matrix as seen in tab. 4.3 (for the G-case). As expected, the constraints in 
the G-case are not as tight as those in the if -case, since we allow the growth 
function to vary in a more general way. 



4.5 Including priors 



As mentioned in section |2.4.1 the magic of Fisher lets us easily include 
priors from other experiments. All we need to do is line up the right rows 
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Figure 4.8 Various parameters are, again, being fixed succes- 
sively (G-case). Here, knowledge about Q m and a% pays off most 
for FOM(iT) (top) and FOM(G) (bottom) respectively. 



and columns of two different Fisher matrices and take their sum, i.e. 

K i-u • i-rr- ( 4 - iQ ) 

We can do this for two measurements of the cosmic microwave background 
radiation. The first one is an already completed experiment, namely WMAP, 
and the second one a planned mission, Planck, both of which use the same 
fiducial model as in this thesis. 



4.5.1 WMAP 



A Gaussian prior, that is a fiducial value Pi with Gaussian error cr(pi), is 
represented simply by a Fisher matrix that has o(pi)~ 2 as its z-th diagonal 
entry (Albrecht et al. : 2006). Gaussian priors can be obtained from recent 



observations. We use the data from the WMAP survey, one of the most 
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a(u m ) 


a{oj b ) 


a(n s ) 


a(Q m ) 




H-case 


0.17 


0.065 


0.049 


0.012 


0.011 


G-case 


0.17 


0.066 


0.059 


0.021 


0.041 



% 


l 


2 


3 


4 


5 


a (hi) 


0.038 


0.015 


0.033 


0.026 


0.053 


a(hi) 


0.044 


0.036 


0.044 


0.049 


0.48 




0.039 


0.023 


0.068 


0.094 


0.88 



Table 4.1 Errors on the cosmological parameters as well as 
hi,hi,gi while fixing r, w , w\, S , 7 and keeping J\f = 5, £ max = 
10 4 , Alg^ = 0.2. 



a(u m ) a{u) h ) a{n s ) a(VL m ) a(w Q ) er(wi) ojj) q-(Eq) er(er 8 ) 
0.23 0.081 0.077 0.045 3.0 9.3 0.67 0.40 0.11 

Table 4.2 Errors on the cosmological parameters while fixing r, 
marginalizing everything else, and keeping A/* = 5, f ma x = 10 4 , 
Alg£ = 0.2 in the fl-case. 



precise measurements of cosmological parameters as of now. The standard 



deviations can be taken from table 3.3 and — while taking into consideration 



the order we are using — the corresponding Fisher matrix thus becomes 

f wmap7 = diag [ 77160 .5 ) 3.42936e6, 5102.04, 

6944.44, 9291.96, 0, 0, 0, 0, 1736.11, 0, . . . , 0] . (4.11) 



4.5.2 Planck 

The Fisher matrix F planck for the planned Planck mission is shown in ta- 



ble 4.4 This matrix needs to be transformed for compatibility with our 



Fisher matrix. The parameter transformation is given by 

0'(0) = (0 5 , 6 6 , const., 6 7 , 1 - 6 3 , 9i, 2 , const., . . . , const) , (4.12) 

where 6' = (u m , ojb, t, n s , Q m , wq, wi, 7, S , as, hi, . . . , gi, . . . ) is the vector in 
parameter space in our formalism, and 6 = (wo, w a , ^de, &k, ^m, k>&, n s ) is 



the vector in parameter space in the formalism of Mukherjee et al. (2008). 
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n b ti 2 


n s 




0"8 


hi 


5.89991e6 


-1.29794e7 


3.04031e6 


1.77766e7 


8.47979e6 


-3.10712e6 


-1.29794e7 


2.85559e7 


-6.68462e6 


-3.91594e7 


-1.86837e7 


6.8399e6 


3.04031e6 


-6.68462e6 


1.57342e6 


9.0384e6 


4.30423e6 


-1.59096e6 


1.77766e7 


-3.91594e7 


9.0384e6 


6.56348e7 


3.13299e7 


-1.06167e7 


8.47979e6 


-1.86837e7 


4.30423e6 


3.13299e7 


1.50269e7 


-5.05321e6 


-3.10712e6 


6.8399e6 


-1.59096e6 


-1.06167e7 


-5.05321e6 


1.99529e6 ... 


-1.94576e6 


4.28365e6 


-9.92991e5 


-7.29185e6 


-3.45161e6 


9.6474e5 


-1.16652e5 


2.56784e5 


-5.94073e4 


-4.72584e5 


-2.22597e5 


4.63944e4 


4.1346e6 


-9.11132e6 


2.09669e6 


1.51661e7 


7.2877e6 


-2.56814e6 


1.91022e6 


-4.20561e6 


9.75038e5 


6.89405e6 


3.2877e6 


-9.21068e5 


4.52229e4 


-9.9446e4 


2.3215e4 


1.77003e5 


8.31533e4 


-1.48989e4 



h 2 

-1.94576e6 
4.28365e6 
-9.92991e5 
-7.29185e6 
-3.45161e6 
9.6474e5 
1.03799e6 
7.90972e4 
-1.55933e6 
-9.6674e5 
-3.22527e4 



h 3 

-1.16652e5 
2.56784e5 
-5.94073e4 
-4.72584e5 
-2.22597e5 
4.63944e4 
7.90972e4 
7.59839e3 
-9.07882e4 
-7.37654e4 
-3.24212e3 



5^ 

4.1346e6 
-9.11132e6 
2.09669e6 
1.51661e7 
7.2877e6 
-2.56814e6 
-1.55933e6 
-9.07882e4 
3.61394e6 
1.48555e6 
3.17958e4 
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1.91022e6 

-4.20561e6 
9.75038e5 
6.89405e6 
3.2877e6 

-9.21068e5 
-9.6674e5 

-7.37654e4 
1.48555e6 
9.18941e5 
3.03478e4 



93_ 

4.52229e4 
-9.9446e4 
2.3215e4 
1.77003e5 
8.31533e4 
-1.48989e4 
-3.22527e4 
-3.24212e3 
3.17958e4 
3.03478e4 
1.4814e3 



Table 4.3 The full Fisher matrix in the G-case for three redshift 
bins. 
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w 


W a 


^DE 




.172276c+06 


.490320e+05 


.674392e+06 


W a 


.490320c+05 


.139551c+05 


.191940e+06 


^DE 


.674392c+06 


.191940e+06 


.263997e+07 


Hfc 


-.208974c+07 


-.594767e+06 


-.818048e+07 . . . 


Mm 


.325219e+07 


.925615e+06 


.127310e+08 


Ub 


-.790504e+07 


-.224987e+07 


-.309450e+08 


n s 


-.549427e+05 


-.156374e+05 


-.215078c+06 





UJm 


^b 


ns 


-.208974e+07 


.325219e+07 


-.790504e+07 


-.549427e+05 


-.594767e+06 


.925615e+06 


-.224987e+07 


-.156374e+05 


-.818048e+07 


.127310e+08 


-.309450e+08 


-.215078c+06 


.253489e+08 


-.394501c+08 


.958892e+08 


.666335e+06 


-,394501e+08 


.633564e+08 


-.147973e+09 


-.501247c+06 


.958892c+08 


-.147973e+09 


.405079e+09 


.219009c+07 


.666335e+06 


-.501247e+06 


.219009e+07 


.242767e+06 



Table 4.4 Fisher matrix F planck for (wq, 



n 



DE- 



k, U} m , Mb, 



derived from the covariance matrix for (R,l a ,Qbh ,n s ) from 



Planck (Mukherjee el al. 2008). 



The Jacobian matrix thus reads 



801 

de 3 



( 



0-1 



1 
10 

\ 1 



1 





.. 







1 


.. 


. 








.. 


. 








.. 


. 








.. 


. 








.. 


. 








.. 


• o/ 









(4.13) 



10+Af(+A0 columns 

The zeros on the right hand side account for the his (or his and g^s), which 
the original Planck Fisher matrix naturally does not contain. The Fisher 
matrix with the PLANCK prior for our use is then 



-Planck 



TcPlanck 



J 7 F 



J. 



(4.14) 



4.5.3 Summary 

We can now examine how including priors from WMAP, Planck, or both 
will affect the constraints on H(z) or G(z). For this, we will choose M = 5 
while marginalizing over all other parameters, although G(z) is fixed when 



calculating the FOM(iJ). The results are summarized in table 4.5 
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WL 


WL+WMAP7 


WL+Planck 


WL+Planck+WMAP7 


FOM(H) 


7636 


13778 


25194 


26096 


FOM(G) 


2865 


5697 


4350 


6805 



Table 4.5 Comparison of the FOM when including different priors 
for the H-case and the G-case, respectively. 



Naturally, the FOM improves in both cases as we include additional pri- 
ors. Unexpectedly, the growth function is more constrained by assuming 
WMAP priors compared to Planck priors. With Planck being one gen- 
eration ahead of WMAP, it should be the other way around. But if we 
take fig. 4^8 into consideration, it is clear that while fixing Q m is most im- 
portant for constraining G, fixing <7g also almost doubles the FOM. Unlike 
the WMAP prior, the Planck prior does not constrain at all, thus even 
though PLANCK constrains fl m much more than WMAP in this parameter 
configuration, it is not as important with regard to constraining the growth 
function. Another observation we can make is that in either case the FOM 
does not increase significantly compared to the next higher value when in- 
cluding the second prior, i.e. the FOM in the if-case, when including WMAP 
and Planck, is almost the same as only including Planck. In the G-case, 
the FOM is almost the same when including both priors as compared to only 
including WMAP. This is to be expected, as both experiments probe the 
cosmic microwave background radiation, except with different sensitivities. 
Therefore, no more cosmological information can be extracted. 



4.6 Comparison to current constraints 

When we try to compare how weak lensing in this context will constrain 
the growth function compared to current constrains, we notice that our con- 



straints are very poor at high redshifts. Rapetti et al. (2010) measured the 
growth index via a combination of X-ray cluster growth data with cluster gas 
mass fraction, type la supernovae, baryon acoustic oscillations and cosmic 
microwave background data, where they found that 7 = 0.551q io- As shown 



in fig. 4.9, these constraints are not very tight when applied to the growth 
function (even when assuming that Q m is known exactly), but the error bars 
from the Fisher matrix analysis hardly seem any better, especially at high 
redshifts. It is interesting, though, to have a model independent confirmation 
of our observations that is comparable with current constraints. 
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0.5 



0.4- 




Figure 4.9 Current constraints on the growth function via the 
growth index by using X-ray cluster growth data. Shown is the 
lcr confidence band of the growth funcion. 



4.7 Possibility of ruling out modified gravity 
models 



According to the data displayed in fig. |4.4[ we get the best results for the 
growth function when using three redshift bins, which is why we shall adopt 
this number for this and the next section. Fig. 4.10| shows the error bars 
on G(z) in this case, along with the growth function as predicted by the 
Starobinsky model with 7s ta = 0.42 (Fu et al, 2010 Starobinsky, 2007) and 



the DGP model with 7dgp by Dvali, Gabadadze and Porrati (Gong, 2008 



Dvali et al. , 2000). We also see the same diagram with priors from WMAP 
and Planck included. From the visual, we recognize a slight disagreement 
at low redshifts while the error bars at high redshifts are too large to make 
any statements about the competing models. 

One caveat: The Starobinsky model actually predicts a redshift depen- 
dent growth index. Trying to fit a power law to the growth rate in this case, 
i.e. 



f = n r 



>Sta 



with a constant 75*0, yields a deviation of up to 10% at 



intermediate redshifts, which indicates that the result cited above is biased. 
A linear expression for the growth index, 7g ta = 70 + Jiz, is more accurate 
but fails at high redshifts with z > 0.5. However, the result, 71 ~ —0.24 
indicates a decreasing growth index (Fu et al, 2010). Thus, in the case of 



the Starobinsky model the power law does not seem to be a good fit and the 
following discussion should be taken with a grain of salt. 
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We now want to attempt to quantify the deviation from these two models. 
Assuming the errors on the growth function are Gaussian, we can compute 
the likelihood function for the growth index 7 as 



£(7) = exp -- ^ — - 2 (4.15) 



using the error estimates from the Fisher matrix and keeping Q m fixed. The 



likelihood is plotted in fig. 4.11 Obviously, the maximum of the likelihood 
is found at the fiducial value 7 = 0.55. In this case, the likelihood can be 
approximated by a Gaussian with good agreement, and we find a standard 
deviation of 00(7) = 0.056 with weak lensing only and 01(7) = 0.037 while 
including priors from WM AP and Planck. Considering that both the DGP 
model as well as the Starobinsky model deviate from the fiducial value by at 
least A7 = 0.13, they could both be ruled out at the 2a confidence level with 
weak lensing only and at the 3a confidence level with priors. This does not 
seem to be entirely conclusive, especially when we consider the uncertainty 
in the 7st a that we mentioned earlier. However, we need to keep in mind that 
these constraints are model independent. 



4.8 Problems and outlook 



It is not surprising that the Fisher matrix formalism can only provide a 
lower bound on the expected error bars. This bound is called the Cramer- 
Rao bound, which expresses the fact that the maximum likelihood estimate 
constrains the model parameters at best with the value given by the Fisher 
matrix (Shao, 2003). We should also keep in mind that when the errors 



become too large, the approximation of the likelihood by a Gaussian may 
lose its validity. 

Besides that, we did not take into account several systematic errors that 
might present itself while observing the weak lensing signal, like the afore- 
mentioned intrinsic alignment. Other potential error sources include: shear 
calibration errors, intrinsic galaxy alignment, and non-Gaussian statistics, 



some of which are discussed in Bernstein (2009) and Albrecht et al. (2009). 



Regarding the source code of the notebook, while it does satisfy all the 
requirements for calculating the Fisher matrix, there certainly is still substan- 
tial room for improvement. One desirable feature would be the abstraction 
of cosmological parameters, such that they would not have to be hard coded 
into the arguments of any function, making it easier to add or remove pa- 
rameters later on. In appendix |AJ large parts of the code were reused, but 
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almost all arguments had to be adjusted manually, a cumbersome and error 
prone job. Another issue, that the runtime on a typical machine is of order 
several hours to a few days with parallelizing part of the code, which could 
perhaps be reduced significantly by examining the step size of interpolated 
functions more carefully. Many of them (and especially the window func- 
tions, a very central part of the computation) only have a single extremum 
and tend to approach zero after that. Limiting the evaluation of sampling 
points specifically to regions where the function takes non-negligible values 
might speed up the computation drastically. From a more structural point 
of view, a more modular setup by using Mathematica packages would be 
desirable to avoid duplicate code fragments, which would make the program 
more expandable, robust, and easier to use. 
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Figure 4.10 The growth function with error bars for the binning 
with the highest FOM. Also shown (from the top): The growth 
function from eq. (3.68) for the DGP model 7dgp = 11/16 ~ 0.69, 
the standard model 7 = 6/11 ~ 0.55 and the Starobinsky model 
with 7sta = 0.42. Top panel: Weak lensing only. Bottom panel: 
Including priors from WMAP and Planck. 
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Figure 4.11 The likelihood for the growth index with a fitted 
Gaussian (dashed line). Indicated are the predicted values for 
the Starobinsky, DGP and standard model. Top panel: Weak 
lensing only. The Gaussian has a standard deviation of a = 0.056. 
Bottom panel: Including priors from WMAP and Planck. Here, 
the standard deviation is a = 0.037. 



Chapter 5 
Conclusion 



We presented and developed the formal basis of the theory of dark energy, 
weak lensing, and statistics in cosmology (in particular the Fisher matrix 
formalism), on which the remainder of this thesis built. Furthermore, we 
described the specifics of the methods of weak lensing tomography using fit- 
ting formulas for the matter power spectrum with non-linear corrections and 
how it is used to compute the observable, the convergence power spectrum. 
It was demonstrated how we avoided assuming an underlying model for the 
Hubble parameter and growth function by using linear interpolation based 
on a choice of redshift bins. We showed how the convergence power spectrum 
is related to the matter power spectrum and how the matter power spectrum 
can be separated into the growth function and the transfer function. We gave 
the fitting formula for the transfer function and how the non-linear correc- 
tions are applied. Finally, we outlined how to implement the computation in 
Mathematica and what kind of challenges might arise in doing so. 

Then we analyzed the behavior of the Fisher matrix when varying differ- 
ent parameters and why it makes sense to use their particular values. We 
showed how to include priors from other experiments and how it improves 
the FOM, a quantity that we defined as a quality measure. We found that 
a small number number of redshift bins is sufficient to extract most of the 
cosmological information from the weak lensing signal - five for best con- 
straints of H(z) and three for G(z). It turned out that the error bars on 
the growth function are most sensitive to the fractional matter density Q m 
and power spectrum amplitude a 8 . The results indicated that the constraints 
that we can expect are comparable to current constraints from several probes 
combined. According to a likelihood analysis, we conclude that future weak 
lensing surveys, using this kind of technique on H(z) and G(z), might be 
able to rule out the Starobinsky and DGP model for modified gravity only 
at a 2(j level when used by itself. Together with priors from surveys of the 



5.0 
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cosmic microwave background like PLANCK, this might improve to a 3a level 
instead. 

It may not have the kind of certainty that meets the ambitious stan- 
dards of physics, but this the price we pay for parameterizing cosmological 
quantities this way. As a model independent approach it represents an im- 
portant check that is not biased towards a model, supplementary to other 
measurements. 

Regardless, the next decade promises to be an exciting one for cosmology 
and dark energy research as new probes and unprecedented high precision 
data will become available from Euclid and similar surveys. Perhaps it is 
weak lensing what will unravel the mystery that is dark energy by providing 
the necessary insights into the evolution of the growth of structure. . . we can 
only imagine what the consequences for physics will be. 
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Appendix A 



Weak lensing Fisher matrix for 
coupled quintessence 



With quintessence being one of the most popular candidates for dark energy 
(Copeland et al, 2006), it is interesting to study the coupling between a 
scalar field and matter. A simple case would manifest in the conservation 
equations with a coupling constant C such that 



rpjjLU 
rpflV 

1 (0);^ 



(A.l) 
(A.2) 



where and Tjf^ describe the energy-momentum tensor of matter and 



the scalar field respectively (see Amendola (1999); Di Porto and Amendola 
( 2008[ )). It is convenient to define a dimensionless coupling constant f3 pro- 
portional to C, i.e. 

^ c 



The evolution of the density contrast 5(a) must then satisfy 



5"(a) + 



m 







(A.3) 



(A.4) 



where a = In a and ' = d/dlna. In an attempt to approximate the solution 
5(a) and considering that G(z) = 5(z), one can try to generalize the solution 
of the original equation (with j3 = 0) phenomenologically as 



G(z) = exp 



fi m (^(l + c/3 2 



dz' 
l + z 



(A.5) 
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p 


a 




h 


n b 


n s 


6.08293e5 


-1.28992e5 


5.48499e6 


3.73233e6 


3.3104e6 


5.02321e5 


-1.28992e5 


2.80677e4 


-1.2008e6 


-8.22476e5 


-7.44588e5 


-1.044e5 


5.48499e6 


-1.2008e6 


5.28593e7 


3.61948e7 


3.27805e7 


4.54416e6 


3.73233e6 


-8.22476e5 


3.61948e7 


2.48333e7 


2.26252e7 


3.06263e6 


3.3104e6 


-7.44588e5 


3.27805e7 


2.26252e7 


2.1013e7 


2.62304e6 


5.02321e5 


-1.044e5 


4.54416e6 


3.06263e6 


2.62304e6 


4.55493e5 



Table A.l The weak lensing Fisher matrix for coupled 
quintessence with M = 5, £ m£C( = 10 , A\g£ = 0.2. 



Di Porto and Amendola (2008) found that a well working least square fit to 



the exact solution, derived numerically, yields the values 



7 = 0.56. 



2.1. 



(A.6) 



Provided that we are given the transfer function for this scenario to plug 



into eq. (3.8 ), we can skip the fitting procedure for the matter power spectrum 



in section 3.2.1 and proceed directly with the non-linear corrections and 
the convergence power spectrum. Now we can run the code with a slightly 
modified growth function and a fiducial model with the parameters 

(f3, a, tt c , h, tt b , n s ) = (0.01, 0.2, 0.226, 0.703, 0.0451, 0.966) , (A.7) 

where a is the exponent in the potential for the scalar field, which we take 
here to be 

V(<f>) = V a <jT a . (A.8) 

Using Af = 5, £ max = 10 4 , A\g£ = 0.2, we find the resulting Fisher to 
be as in tab. |A.l| We can now show the correlation of all possible parameter 
pairs by marginalizing over all other parameters such that the correlation 
matrix takes the form 

C = ( ° l pa f 2 ) . (A.9) 

Then we can plot ellipses whose semiaxes are proportional to the square root 
of the eigenvalues of C, i.e. 

<2 = I (oi + 4 T sJat-2alal + Ap^alai + at\ (A.10) 
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Weak lensing Fisher matrix for coupled quintessence 



with an angle with respect to the coordinate axes of 



tan 2a 



2 pa i cr 2 
a\-a\ 



(A.ll) 



The semiaxes need to be scaled by a factor of 1.51, 2.49, 3.44 so that the area 
of the ellipsis corresponds to the la, 2a, 3a probability content respectively 



(A&T, ch. 13.3). The result can be seen in fig. A.l and A. 2 (Amendola et al. 



2011). These plots enable us to judge how different surveys can complement 
each other, e.g. we could overlay ellipses derived from CMB (if we had it 
available) data and see how much the confidence regions overlap. Ideally, we 
would want to see thin, perpendicular ellipses. 
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Figure A.l Confidence ellipses for dark energy parameters in the 
coupled quintessence model. Shown are the ellipses for the 68%, 
95% and 99% confidence level. 




Figure A. 2 Confidence ellipses for dark energy parameters in the 
coupled quintessence model (continued). 
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